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We report on a set of simulations performed at several values of the lattice spacing on the Cray T3D at 
Edinburgh. Different methods to extract the lattice scale from the static quark potential are discussed. 



Over the past three years, the UKQCD Collab- 
oration has produced a great many of results in 
lattice QCD using data that were generated on 
the 64-node, i860-based Meiko Computing Sur- 
face ( "Maxwell" ) . Now the collaboration has em- 
barked on a series of new simulations that are 
carried out on the recently installed Cray T3D at 
Edinburgh. This machine is based on the MIMD 
concept, and currently its version at Edinburgh 
consists of 128 nodes, each containing two DEC- 
Alpha processors capable of a peak speed of 300 
MFlops and with 128 MBytes of memory. Thus, 
the peak speed of the whole system is 38 GFlops 
with a total memory of 16 GBytes. The nodes 
form a Torus in 3 Dimensions with the possibility 
of bidirectional communication. A fairly unusual 
feature is the powerful front-end, a Cray Y-MP 
of which there is a two-processor version at Ed- 
inburgh. The Y-MP is Hnked to 200 GBytes of 
disk space as well as to a mass data store with a 
capacity of 1 TByte which is to be upgraded to 
about 20 TBytes in the future. 

In order to ensure portability as well as a long 
lifetime of the code used in our simulations of 
quenched lattice QCD, UKQCD has decided to 
employ standard Fortran?? with PVM for mes- 
sage passing as our programming model [1]. In 
addition, a High-Performance Fortran version has 
been developed, which however has not been used 
in any of the production runs so far. 

The current sustained speed is about 4 GFlops 
for the pure gauge code and about ? GFlops for 
the propagator code (measured on the full 256 



Table 1 

Lattices simulated, showing the ratio of overre- 
laxed to Cabibbo-Marinari sweeps, the number 
of separating sweeps and the number of configu- 
rations analysed for this study. 
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processor partition). The average I/O speed is 
20 MBytes/sec. It should be emphasised that the 
code is still being tuned towards better perfor- 
mance concerning both speed and I/O. 

The availability of a Y-MP as front-end makes 
it possible to perform some part of the analysis 
code on the vector processors which takes some 
load off the main machine. This idea is exploited 
in the computation of gluonic observables. 

The runs on the Cray T3D represent a contin- 
uation of UKQCD's earlier simulations, now ex- 
tended to a range of /3- values and high statistics. 
A summary of lattice sizes and /3-values is pre- 
sented in Tab. 1. Our ultimate goal is to study 
hadron phenomenology including a quantifica- 
tion of residual lattice artefacts using the 0{a)- 
improved Sheikholeslami-Wohlert action [2,3]. 

The gauge configurations were generated on 
the T3D using the Cabibbo-Marinari heatbath 
method plus a number of overrelaxation sweeps. 
They were subsequently analysed on the Y-MP 
front-end to measure the static quark potential, 
Polyakov loops and the two lowest glueball states. 
These measurements mainly serve to set the lat- 
tice scale, using either the string tension ^/K or 
the hadronic scale ro whose merits are discussed 



in [4]. 

The static quark potential V{r) was extracted 
by measuring "fuzzed" Wilson loops W{t,r), 
following the procedure outlined in [4]. At 
/3 = 6.0, 6.2 the orientations (1,0,0), (1,1,0) and 
(1,1,1) were computed, whereas at /? = 6.4 only 
(1,0,0) was considered. The fuzzed links were con- 
structed with link/staple ratio of 2:1 and with 
L — 2 and L + 8 iterations of the fuzzing algo- 
rithm (L being the spatial lattice size). The two 
fuzzing levels were taken as a 2 x 2 variational ba- 
sis to extract the eigenvalues and -vectors of the 
generalised eigenvalue equation 



W,,it,r)v{r)f^ = \k{t,r)W,,it,r)v{r)f\ (1) 



The eigenvector w(r)(^) corresponding to \i{to,r) 
aX to = 3 was used to project onto the approxi- 
mate ground state. The resulting correlator was 
then fitted to both single and double exponentials 
for timeslices up to i = 8. Unless stated other- 
wise the following analysis is based on the results 
from double exponential fits. Statistical errors 
are obtained from a bootstrap analysis based on 
100 bootstrap samples. 

In order to obtain the string tension vK, the 
data for the potential were fitted according to 



V{r)=Vo + K\r\-eGL(r), 



(2) 



where Giir) is the lattice Greens function for 
one-gluon exchange. Our best estimates for the 
fitted parameters are shown in Tab. 2. We have 
included an estimate of systematic errors based 
on the spread obtained from varying the ground 
state potential and the minimum distance in- 
cluded in the fit to eq. (2). The value for the 
string tension are in good agreement with those 
quoted in refs. [5,6]. In the last column we quote 
the values obtained from determining the string 
tension by measuring the correlation of Polyakov 
loops (i.e. the torelon mass). The signal for the 
torelon is still relatively poor, but the values are 
in good agreement with the Wilson loop data and 
also with ref . [7] . 

As discussed in [4], an alternative lattice scale 
is the distance ro at which 



This relation corresponds to ro — 0.5 fm in phe- 
nomenological non-relativistic potential models 
describing quarkonia. In order to compute ro , we 
calculated the tree- level improved force F^rj) [4] 
according to 

F,<ri) = \d\-^Vir)-V{r-d)), (4) 

r, = \d\(GL{r)-GL{r-d)y\ (5) 

ro was then obtained by linearly interpolating the 
force around F{ri) rj = 1.65 [4]. Tab. 2 shows the 
results including the estimates of systematic er- 
rors from the spread of values obtained using dif- 
ferent points in the interpolation and by varying 
the estimate for the ground state potential. 

Our results for ro and ro vK are consistent 
with those quoted in [9] and [8]. However, we 
note that our statistical errors, in particular at 
/? = 6.4, are smaller. 

Using the values for ro we can now scale the re- 
sults at all three /3- values and convert into phys- 
ical units using ro = 0.5 fm. This is shown in 
Fig. 1. The data can be compared to the predic- 
tion from the string formula 



F{r)rl =1.65+ -(-^-1 
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(6) 



F{ro)rl = 1.65. 



(3) 



which is derived from (2) using the constraint (3) 
and e = 7r/12. This gives a remarkably good de- 
scription of the data (solid line) . The linearity of 
the force around 0.5 fm illustrates the fact that ro 
can be obtained reliably since one is not required 
to know the detailed behaviour of the curve at 
very small or large \r\. 

The /3-dependence of the results for ro can be 
compared to perturbation theory using the two- 
loop /3-function for the gauge coupling. The value 
of ro at /3 = 6.4 was used as input, and as was 
taken from ref. [10]. Fig. 2 compares the two- 
loop prediction with the one-loop formula. It is 
seen that the two-loop formula describes the data 
rather well, strengthening our confidence in the 
consistency of our results. 

Using ro = 0.5 fm and Vk = 440 MeV we can 
now give estimates for a~^ GeV. The results in 
Tab. 3 again include our estimates of systematic 
uncertanties. 



Table 2 

The fitted parameters for the potential, eq. (2), results for rg, and the string tension from the Polyakov 

loop ("PL"). The first error is statistical, the second is an estimate of systematic uncertainties. 
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Figure 1. The force at all /3- values, obtained from 
scaling the results using tq. Statistical errors on 
ro are included. 



Table 3 

Estimates for a~^ [GeV] from ro and the string 
tension obtained from the potential ("V") and 
the Polyakov loop ("PL"). 
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Our results represent a first step in our at- 
tempts to perform a controlled extrapolation to 
the continuum limit. The analysis of the static 
quark potential shows a consistent picture for the 
three values of /? considered. In particular, it is 
confirmed that ro can be obtained reliably from 
the linear part of the interquark force. These find- 
ings need to be corroborated with higher statis- 
tics. 

I wish to thank S. Booth, D. Henty, M. Liischer, 
C. Michael, D. Richards, R. Sommer and M. Teper 
for useful discussions and suggestions. 
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Figure 2. The values for ro at different /? com- 
pared to the predictions of 2-loop (solid line) and 
1-loop (dashed line) perturbation theory. 
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